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Abstract 

We report the discovery of a fifth image in the large separation lensed quasar system SDSS J1004+4112. 
A faint point source located 0"2 from the center of the brightest galaxy in the lensing cluster is detected 
in images taken with the Advanced Camera for Surveys (ACS) and the Near Infrared Camera and Multi- 
Object Spectrometer (NICMOS) on the Hubble Space Telescope. The flux ratio between the point source 
and the brightest lensed component in the ACS image is similar to that in the NICMOS image. The 
location and brightness of the point source are consistent with lens model predictions for a lensed image. 
We therefore conclude that the point source is likely to be a fifth image of the source quasar. In addition, 
the NICMOS image reveals the lensed host galaxy of the source quasar, which can strongly constrain the 
structure of the lensing critical curves and thereby the mass distribution of the lensing cluster. 
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1. Introduction 

Gravitational lensing is a unique tool for exploring the 
distribution of matter, particularly that of dark matter. 
The recently discovered largest separation lensed quasar, 
SDSS J1004+4112 (Inada et al. 2003; Oguri et al. 2004), 
has opened a new window for probing dark matter dis¬ 
tributions in the universe. The quasar was discovered in 
the Sloan Digital Sky Survey (SDSS; York et al. 2000; 
Abazajian et al. 2004), and the lensing hypothesis was 
confirmed by subsequent observations with the Subaru 
8.2-m telescope and the Keck telescope. The system con¬ 
sists of four lensed quasar components {%' =18.5, 18.9, 
19.4, and 20.1) at 2 = 1.734 and the maximum separa¬ 
tion angle between the lensed images is 14"62. The lens¬ 
ing object must be a massive object, such as a cluster of 

* Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. These 
observations are associated with HST program 9744. 


galaxies, to produce such a large image separation; indeed, 
we have identified a. z = 0.68 cluster centered among the 
four lensed images. The discovery of a single, cluster-size 
lensed quasar among the current SDSS quasars is consis¬ 
tent with the theoretical expectation of lensing based on 
the cold dark matter model (Oguri et al. 2004; Oguri & 
Keeton 2004). 

SDSS J1004+4112 is unique in the sense that 1) the 
quadruple images place robust constraints on the inner¬ 
most region of the lensing cluster, and 2) the lensing clus¬ 
ter is a strong lensing selected cluster of galaxies. These 
features indicate that the mass modeling of this lens sys¬ 
tem may offer valuable information on the structure of 
clusters of galaxies. The first attempt at modeling the 
system with various parametric models revealed the elon¬ 
gated and complicated mass distribution of the lensing 
cluster (Oguri et al. 2004). Williams & Saha (2004) re¬ 
cently studied this system using a free-form reconstruc¬ 
tion technique and reached similar conclusions. However, 
important degeneracies between different models remain, 
and further follow-up observations are required to deter- 
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mine the mass distribution more precisely. 

Observations using the Hubble Space Telescope ( HST) 
can offer such new data. In particular, high-resolution 
HST images could be quite effective at detecting lensed 
images of the quasar host galaxy, arc or arclet images 
of other lensed sources, and perhaps a central or “odd” 
image of the lensed quasar system that would be ex¬ 
pected theoretically (e.g., Burke 1981; Rusin 2002). A 
central image is especially useful for providing tight con¬ 
straints on the central mass distribution of the lensing 
object. This was indeed demonstrated by the first discov¬ 
ery of a central image in a lensed quasar system; Winn, 
Rusin, & Kochanek (2004) showed that the central image 
of PMN J1632—0033 requires /? = 1.91 ± 0.02 (2a confi¬ 
dence) when a power-law density profile p(r) oc r“' 9 is as¬ 
sumed. Possible central images have also been identified in 
lensed arc systems, such as CL 0024+1654 (Colley, Tyson, 
& Turner 1996), MS 2137.3-2353 (Gavazzi et al. 2003), 
and A1689 (Broadhurst et al. 2005), and they also provide 
important constraints on mass models (see, e.g., Gavazzi 
et al. 2003). 

In this letter , we present an identification of a 
fifth image of the lensed quasar with the Advanced 
Camera for Surveys (ACS; Clampin 2000) and the 
Near Infrared Camera and Multi-Object Spectrometer 
(NICMOS; Thompson 1992) installed on the HST. In ad¬ 
dition, we report unambiguous detection of the lensed host 
galaxy in the NICMOS image. 

2. The ACS Observation 

An ACS observation (with the Wide Field Channel) in 
the F814W filter («/-band) was conducted on 2004 April 
28, under the program “ HST Imaging of Gravitational 
Lenses” (GO-9744, PI C. Kochanek). The observation 
consisted of five dithered exposures taken in ACCUM 
mode. The total exposure time was 405 seconds. The 
reduced (drizzled and calibrated) images were extracted 
using the CALACS pipeline (Hack 1999), which includes 
the PyDrizzle algorithm. We further rejected cosmic rays 
using the L.A.Cosmic package (van Dokkum 2001) in the 
drizzled image. A 40" x 40" subsection of a median of 
the five dithered images is shown in Figure 1. Following 
Figure 9 of Oguri et al. (2004), the four lensed components 
are denoted as “A-D”, and three central bright galaxies 
of the lensing cluster are denoted as “G1-G3”. The red- 
shifts of galaxies G1-G3 are 0.680, 0.675, and 0.675, re¬ 
spectively (Oguri et al. 2004). The relative positions of 
components A D were calculated by a single Gaussian fit, 
and the flux ratios of components A D were estimated by 
fitting the PSF stars produced by the Tiny Tim software 
(version 6.1a; Krist & Hook 2003), in the drizzled (and 
cosmic ray rejected) image. The PSF of a quasar was 
constructed with an a v = —0.5 power law spectrum in the 
F814W wavelength region (corresponding to ~3000 A in 
the rest frame). The relative positions derived from the 
ACS (F814W) are consistent with those derived from the 
Subaru i' band image (Oguri et al. 2004) within ~5cr. The 
four components are all unsaturated, and the AB magni- 


Table 1. Relative Positions and Flux Ratios of 
SDSS J11004+4112 in the HST ACS Image 


Object 

x[arcsec]* 

y[arcsec]* 

Flux Ratio 1, 

A 

0.0000+0.001 

0.0000+0.001 

1.000 

B 

-1.317+0.002 

3.532+0.002 

0.732 

C 

11.039+0.002 

-4.492+0.002 

0.346 

D 

8.399+0.004 

9.707+0.004 

0.207 

E 

7.197+0.009 

4.603+0.009 

0.003 

Gl 

7.114+0.030 

4.409+0.030 



respectively. 

t Errors of fitting by quasar PSFs are about 10%. The error of com¬ 
ponent E might be much larger due to the over-subtraction of Gl. 


tude of component A in the F814W filter is estimated to 
be 18.4. The position of component Gl was extracted by 
the Source Extractor algorithm (Bertin & Arnouts 1996). 
The results are summarized in Table 1. 

3. The NICMOS Observation 

The NICMOS imaging observation was also conducted 
under the same HST program on 2004 October 9. The 
observation consists of four dithered exposures taken in 
MULTIACCUM mode, using the F160W filter («iF- 
band). The exposure time was 640 sec for two of the 
exposures and 704 sec for the other two. The calibration 
was extracted by the CALNICA pipeline, and the cen¬ 
tral bad columns of each dithered image were corrected 
by linear interpolation. The combined image is shown in 
Figure 2. First, we confirm a galaxy near component A 
(marked as G4 in Figure 2), which may host the star or 
stars responsible for microlensing of the broad emission 
line region (Richards et al. 2004). In addition, extended 
emission is clearly seen around components A, B and C. 
Although such extensions are also seen in the ACS im¬ 
ages (see Figure 1), their existence is much more robust 
in the NICMOS image. The fact that these extensions are 
obvious in the F160W (near-infrared) image and faint in 
the F814W (optical) image, and that the distortions agree 
with the theoretical critical curves (see Figure 17 of Oguri 
et al. 2004), demonstrates that the extended flux is due to 
the lensed host galaxy of the source quasar. These images 
provide many new constraints on lens models that will sig¬ 
nificantly improve our ability to determine the mass dis¬ 
tribution of the lensing cluster. They do make lens mod¬ 
eling more computationally intensive (because one must 
account for the intrinsic shape of the host galaxy, and also 
for the effects of the point spread function), so we defer 
detailed modeling of the arcs to a subsequent paper. 

4. Fifth Image 

Of interest is the existence of a point source near the 
center of Gl. The left panel in Figure 3 shows the ACS im¬ 
age of galaxy Gl; what appears to be an unresolved source 
is clearly seen approximately 0"2 northwest of the center 
of Gl. This feature is neither a bad pixel nor a cosmic ray; 
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Fig. 1. 40” X 40” subsection of the median combined ACS image of SDSS J1004+4112. The four quasar images are labeled A-D. The 
galaxies of the lensing cluster labeled G1-G3 have redshifts of 0.680, 0.675, and 0.675, respectively. The pixel scale is approximately 
0”05 pixel” 1 . 


the source is seen in all dithered images. The right panel 
in Figure 3 displays the image after subtracting the signal 
from G1 (modeled with the GALFIT package of Peng et 
al. 2002). Due to the existence of the unresolved source 
near the center of G1 (the peak flux of this unresolved 
source is almost same as that of Gl), G1 was slightly 
over-subtracted. However, we can see a new source, la¬ 
beled E, in the subtracted image. This object is classified 
as a point source by the Source Extractor algorithm. The 
position and brightness of E, based on a single Gaussian 
fit to the data, are given in Table 1. 

We find component E also in the NICMOS image; the 
left panel in Figure 4 is the NICMOS image of galaxy 
Gl. We subtracted the signal from Gl using the GALFIT 
package, which is shown in the right panel of Figure 4. 
We confirm component E in the subtracted image, al¬ 
though Gl was slightly over-subtracted as in the ACS 
image. Measuring the flux of component E with a sin¬ 
gle Gaussian fit, we find that the flux ratios between E 
and A (E/A) in the ACS and NICMOS images are 0.003 
and 0.004, respectively. This remarkable agreement of the 
flux ratios supports the idea that the point source is a fifth 
image of the lensed quasar. 

To test this hypothesis, we have refined the lens models 
presented in Oguri et al. (2004) using the more precise 
HST data. The models consist of a singular isothermal 
ellipsoid mass distribution for galaxy Gl, and an NFW 
(Navarro, Frenk, & White 1997) elliptical potential for the 
cluster. We demand that the models reproduce the rela¬ 


tive positions and brightnesses of components A-D (and 
the relative position of Gl), as well as the position of 
component E; we do not use the brightness of E as a con¬ 
straint, because we want to see what the models predict. 
Adopting the same approach as Oguri et al. (2004), we 
use the lensmodel software (Keeton 2001) for Monte Carlo 
sampling of the parameter space. There is a wide range 
of models consistent with the data, indicating that it is 
not difficult to produce a 5-image lens matching the con¬ 
figuration of components A -E, and that significant model 
degeneracies remain. More interesting are the model pre¬ 
dictions for the flux ratio between the fifth image and 
component A, shown in Figure 5. The predictions span 
a remarkable nine orders of magnitude, from E/A ^0.1 
down to E/A ~ 10 -10 , but a significant fraction predict 
E/A in the range 0.001-0.01, consistent with the observed 
value. In other words, there are (many) reasonable mod¬ 
els that can fit all of the HST data under the hypothesis 
that E is a fifth image; conversely, the observed properties 
of E are highly compatible with that hypothesis. 

Given the enormous range of model predictions for the 
brightness of E, it appears that the observed brightness 
offers strong constraints on the models. We caution that 
one must be careful in using the brightness of E as a con¬ 
straint, because its measured flux could be contaminated 
by improper subtraction of the galaxy or by physical ef¬ 
fects such as microlensing or extinction. Nevertheless, the 
range of predictions is so large that even conservative es¬ 
timates of systematic uncertainties should still yield very 
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Fig. 2. The combined NICMOS image of SDSS J1004+4112. 
The pixel scale is approximately 0^075 pixel -1 . The lensed 
host galaxy is seen more prominently in this NICMOS image 
than in the ACS image. In addition, we can see G4 near 
component A. 


Fig. 3. Left: The ACS image of central galaxy Gl. Right: 
The ACS image after subtracting galaxy Gl. The residual 
image clearly shows a stellar object near the center of Gl. 

interesting results. As an example, to test the ability 
of E to constrain the central density profile of Gl, we 
switch from an isothermal model to a more general power 
law density profile p{r) oc r~^ for the galaxy (for com¬ 
putational simplicity, we now assume that the potential , 
rather than the density, has elliptical symmetry). We find 
that there is a wide range of models with 1.6 < (3 < 2.0 
that predict 0.001 < E/A < 0.01. By contrast, all of the 
models we examined with (3 > 2.1 predict E/A < 10 -10 , 
which grossly contradicts even a conservative reading of 
the data. In other words, the observed properties of E im¬ 
ply that the galaxy mass distribution cannot be steeper 
than isothermal (i.e., (3 < 2). This upper bound is similar 
to that found by Winn, Rusin, & Kochanek (2004) from 
the central image in PMN J1632—0033 (specifically, they 
found (3= 1.91 ±0.02). However, SDSS J1004+4112 dif¬ 
fers from PMN J1632—0033 in that we do not obtain any 



Fig. 4. Left: The NICMOS image of central galaxy Gl. 
Right: The NICMOS image after subtracting galaxy Gl. The 
stellar object near the center of Gl can be also seen in the 
NICMOS image. 

lower bound on /3, at least over the range 1.6 < f3 < 2.0 
that we have explored so far. Apparently the complexity 
of the SDSS J1004+4112 lens potential, with a cluster in 
addition to the galaxy, prevents a unique measurement 
of the value of (3. Nevertheless, it is clear that compo¬ 
nent E provides important new constraints on the mass 
distribution of this interesting lens system. 

5. Summary 

We have presented the HST ACS and NICMOS images 
of SDSS J1004+4112, which reveal a fifth image of the 
lensed quasar core. The fifth image offers a unique probe 
of the mass distribution of the cluster core. Deep spectro¬ 
scopic observations of component E with large telescopes, 
such as the Subaru Telescope, offer the best prospect for 
the final confirmation that it is a lensed quasar image. In 
the NICMOS image, we also found unambiguous evidence 
of the lensed host galaxy of the source quasar. These 
extended images provide strong additional constraints on 
mass models of the lensing cluster, which are expected to 
break degeneracies seen in the modeling studies to date. A 
detailed analysis of lens models including the host galaxy 
images is underway and will be presented elsewhere. 
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Fig. 5. Frequency (fraction of models) of the predicted flux 
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observed positions and brightnesses of components A-D (and 
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ACS) with 50% error. 

Mexico State University, University of Pittsburgh, 
University of Portsmouth, Princeton University, the 
United States Naval Observatory, and the University of 
Washington. 

References 

Abazajian, K. 2004, AJ, 128, 502 
Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 
Broadhurst, T., et al. 2005, ApJ, in press (astro-ph/0409132) 
Burke, W. L. 1981, ApJL, 244, LI 
Clampin, M. et al. 2000, SPIE, 401, 344 
Colley, W. N., Tyson, J. A., & Turner, E. L. 1996, ApJL, 461, 
L83 

van Dokkum P. G. 2001, PASP, 113, 1420 
Gavazzi, R., Fort, B., Mellier, Y., Pello, R., & Dantel-Fort, M. 
2003, A&A, 403, 11 

Hack, W. J. 1999, CALACS Operation and Implementation, 
ISR ACS-99-03 

Inada, N., et al. 2003, Nature, 426, 810 
Keeton, C. R. 2001, preprint (astro-ph/0102340) 

Krist, J. E. & Hook, R. N. 2003, The Tiny Tim User’s Guide, 
Version 6.1a 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 
490, 493 

Oguri, M., et al. 2004, ApJ, 605, 78 
Oguri, M., & Keeton, C. R. 2004, ApJ, 610, 663 
Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2002, AJ, 
124, 266 

Richards, G. T., et al. 2004, ApJ, 610, 679 
























